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The Einstein-Straus spacetime describes a nonrotating black hole immersed in a matter-dominated 
cosmology. It is constructed by scooping out a spherical ball of the dust and replacing it with a 
vacuum region containing a black hole of the same mass. The metric is smooth at the boundary, 
which is comoving with the rest of the universe. We study the evolution of a massless scalar field 
in the Einstein-Straus spacetime, with a special emphasis on its late-time behavior. This is done 
by numerically integrating the scalar wave equation in a double-null coordinate system that covers 
both portions (vacuum and dust) of the spacetime. We show that the field's evolution is governed 
mostly by the strong concentration of curvature near the black hole, and the discontinuity in the 
dust's mass density at the boundary; these give rise to a rather complex behavior at late times. 
Contrary to what it would do in an asymptotically-flat spacetime, the field does not decay in time 
according to an inverse power-law. 




I. INTRODUCTION 

The dynamics of radiative fields in spacetimes contain- 
ing a black hole has been the subject of many investi- 
gations, some analytical some numerical 
some restricted to nonrotating black holes pJ 
and some devoted to rotating black holes 
Most of these studies were concerned with black holes 
immersed in an asymptotically-flat spacetime. In such 
cases, radiative dynamics always proceeds in the same 
three stages. First, an outburst of radiation is emitted, 
as determined by the initial conditions imposed on the 
radiative field; most of this radiation proceeds directly 
to infinity, with only mild distortions produced by the 
spacetime curvature. Second, the radiation begins to os- 
cillate, with frequencies and damping times determined 
by the properties of the black hole; these oscillations are 
produced by the part of the initial outburst that does 
not go directly to infinity, but interacts with the high 
concentration of curvature near the black hole. Third, 
the oscillations stop, and the field decays monotonically, 
as an inverse power-law in time; the power is determined 
by the radiation's multipole order I. 

A number of authors have studied radiative dynam- 
ics in black-hole spacetimes that are not asymptotically 
flat. Brady, Chambers, Krivan, and Laguna |lq| consid- 
ered the evolution of a massless scalar field in Reissner- 
Nordstom-de Sitter spacetime, which describes a charged 
black hole immersed in an exponentially expanding uni- 
verse. They found that while the first two stages of ra- 
diative dynamics are not affected by the different con- 
ditions at infinity, the third one is: at late times the 
field decays exponentially, not as an inverse power-law. 
This new type of behavior was explored more thoroughly 
by Brady, Chambers, Laarakkers, and Poisson J!?]], who 
showed that the decay constant can be complex if the 
scalar field is nonminimally coupled to the curvature: the 



field oscillates as it decays exponentially. While these au- 
thors were concerned with black-hole spacetimes with a 
positive cosmological constant, other authors have exam- 
ined radiative dynamics in the presence of a negative cos- 
mological constant p8|-pT| ; they also find an exponential 
decay. 

The literature reviewed in the preceding paragraph 
indicates that inverse power-law decay is a property of 
asymptotically- flat spacetimes, and that different behav- 
iors should be expected in black-hole spacetimes that 
have different asymptotic properties. This conclusion 
has been reinforced by a number of analytical studies 
H||. Our goal with this paper is to contribute fur- 
ther to the exploration of radiative dynamics in black- 
hole spacetimes that are not asymptotically flat. Here 
we consider the evolution of a massless scalar field in 
the spacetime of a nonrotating black hole immersed in a 
matter-dominated cosmology. We investigate the effect 
of the cosmological asymptotics on the late-time behavior 
of the radiation. 

To construct the spacetime we start with a homoge- 
neous, isotropic, dust-filled universe and we remove from 
it a spherical ball of the dust, which we replace by a vac- 
uum region containing a black hole of the same mass. 
The boundary between the vacuum and the dust is co- 
moving with the rest of the universe. Its dynamics is 
determined by the curvature of the spatial sections: it 
expands forever for zero and negative curvature, and it 
eventually recontracts for positive curvature. This con- 
struction was first carried out by Einstein and Straus 
in 1945 p2fl . We reproduce it in Sec. II, where we also 
introduce a single coordinate system to cover both re- 
gions (vacuum and dust) of the spacetime. In Sec. Ill we 
write down the second-order partial differential equation 
that governs the evolution of a massless scalar field <& 
in this spacetime. After separation of the variables, the 
field equation reduces to a two-dimensional wave equa- 
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tion with an effective potential. This wave equation is 
numerically integrated in Sec. IV, where we present our 
main results. 

Radiative dynamics in the Einstein-Straus spacetimc 
is governed mostly by the strong concentration of cur- 
vature near the black hole, and the discontinuity in the 
dust's mass density at the boundary. These features are 
encoded in the effective potential, which possesses a fairly 
well localized barrier near the black hole's event horizon, 
and a jump discontinuity at the boundary. A wave prop- 
agating in the presence of this potential will be partially 
transmitted and reflected each time it encounters the bar- 
rier or the jump. These processes give rise to a rather 
complex late-time behavior for the radiation — the field 
does not decay according to an inverse power-law. We 
describe this behavior in detail in Sec. IV. In Sec. V we 
seek insight into the multiple transmissions and reflec- 
tions that take place between the event horizon and the 
boundary by working through a simple toy model. 

We believe that radiative dynamics in the Einstein- 
Straus spacetime provides an excellent illustration of 
the general fact that inverse power-law falloff can occur 
only in asymptotically-flat spacetimes. This article con- 
tributes to an ongoing effort to explore the rich spectrum 
of behaviors that results when asymptotic flatness is re- 
placed by different asymptotic conditions. 

Throughout the paper we work in relativistic units, in 
which c = G = 1 . 



II. EINSTEIN-STRAUS SPACETIME 

The Einstein-Straus spacetime p2| represents a non- 
rotating black hole immersed in a cosmological universe 
containing a pressureless fluid (which we shall call dust). 
It is constructed by scooping out a spherical ball of the 
dust and replacing it with a vacuum region containing 
a black hole of the same mass. The mathematical de- 
scription of the spacetime involves two metrics joined at 
a common boundary E. The metric inside the boundary 
is given by the Schwarzschild solution, the metric out- 
side is one of the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metrics, and £ is a three-dimensional hypersur- 
face generated by timelikc geodesies of both spacetimes. 

The boundary is comoving with the universe; it ex- 
pands forever if the universe has spatial sections that are 
flat (k = 0) or have negative curvature (k = —1), or 
it expands to a maximum radius and recontracts if the 
spatial sections have positive curvature (k = 1). As seen 
from the Schwarzschild side, the boundary is either grav- 
itationally bound to the black hole (if k = 1), or it is 
unbound or just marginally bound (if k = — 1 or k = 0, 
respectively). The spacetime is depicted, by means of a 
conformal diagram, in Fig. 1; in this figure the boundary 
is assumed to expand forever. 



black-hole singularity 




past singularity 



FIG. 1. The Einstein- Straus spacetime, as represented by 
a conformal diagram. The Schwarzschild region is on the 
left-hand side of the boundary, and the FLRW region in on 
the right-hand side. It is assumed that the cosmological region 
possesses flat spatial sections, so that the boundary expands 
forever. The portion of the spacetime bounded by the past 
singularity, the future horizon, and future null infinity can be 
covered by the u and v coordinates introduced in the text. In 
the diagram, u increases at 45 degrees toward the left, while v 
increases at 45 degrees toward the right. The future horizon is 
located at u — oo, future null infinity at v = oo, the boundary 
at v — u — 2xo, and the past horizon at v = 3\o- 

A. Joining the metrics at the boundary 

The metric inside S is expressed in the usual 
Schwarzschild coordinates as 

dsf n = -fdt 2 + f- 1 dr 2 +r 2 dfl 2 , (2.1) 

where / = 1 - 2M/r and dQ 2 = d0 2 + sin 2 Ode/) 2 ; M is 
the grav itational mass of the black hole. The metric of 
Eq. ( |2.l| ) is valid for r < ts{t), where is the changing 
radius of the timelike hypersurface E; r is proper time 
for observers comoving with the boundary. Because the 
hypersurface is generated by timelike geodesies, its equa- 
tions of motion are t = E/ f and r 2 + f = E 2 . Here, an 
overdot indicates differentiation with respect to r, and E 
is the usual conserved-energy parameter. 

For bound motion (E < 1), the solution to the equa- 
tions of motion can be expressed in the parametric form 

rs(??) = -r max (l - cos 77), 

(2.2) 

r (v) = 2 r maxV r max/2Af(?7 - SUIT?), 

in terms of a conformal-time parameter 77; the constant 
r max = 2M/ (1 — E 2 ) represents the boundary's maximum 
radius. For unbound motion (E > 1), the solution is 



2 



rv(v) = 2 r o(coshr/ - 1), 



(2.3) 



T (v) = ^ r ° V r o/2M(sinh?7 - jj), 



where r$ = 2M/(E 2 — 1). Finally, for marginally bound 
motion [E = 1) we have 



_ 9M f 2t 
~ ~2~ I 9M 
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The metric outside E is expressed as 



ds 2 



a 2 (77)[-dr/ 2 +d X 2 + s 2 (x) dfi 2 



(2.4) 



(2.5) 



where 77 is conformal time, x a radial coordinate, a(rj) 
the cosmological scale factor, and s(x) a function that 
determines the geometry of the spatial sections: 



sinx k = 1 

X *: = 

sinhx k = — 1 



(2.6) 



The metric of Eq. (2.5) is valid for x > Xo, with xo de- 
noting the comoving position of E. The behavior of the 
scale factor is determined by the spatial curvature and 
the matter content of the universe. Here we assume that 
the universe is matter dominated, with a stress-energy 
tensor given by T a @ = pu a vfi , where p is the dust's mass 
density and u a = a~ 1 dx a /drj its four- velocity. Conser- 
vation of energy implies that pa 3 — constant = 3C/(8tt). 
By virtue of the Einstein field equations, the scale factor 
must satisfy a' 2 + ka 2 = Ca, in which a prime indicates 
differentiation with respect to 77. Proper time for comov- 
ing observers is obtained by integrating dr — a(rj) dr\. 

For spatial sections with positive curvature [k = 1), 
the scale factor and the proper time are given by 

C C 
a (v) = -w (1 ~ cos?7), t(t)) = — (77- sin 77), (2.7) 

where C — 87rpa 3 /3 is a constant. For spatial sections 
with negative curvature (k = —1), we have instead 

C C 
a(v) = — (cosh?7- 1), t(?7) = — (sinh 77 - 77). (2.8) 

Finally, for flat spatial sections (k = 0) we have 



(2.9) 



The metrics of Eq. (2.1) and (2.2) must be matched at 
their common boundary p3[ , the timelike hypersurface 
S described by r = rs(r) on the Schwarzschild side, and 
X = Xo 011 the FLRW side. The boundary's induced 
metric, as seen from the interior, is 



[r s (r)] 2 dQ 2 



(2.10) 



As seen from the exterior, it is 

dsl = -dr 2 + [a(r)s( Xo )] 2 dQ 2 



(2.11) 



In both cases the induced metric is expressed in terms 
of the hypersurface's intrinsic coordinates (r, 9, </>), with 
t denoting proper time for observers comoving with the 
hypersurface. Equality of the induced metric on both 
sides of the hypersurface implies 



rs(r) = o(t)s(xo), 



(2.12) 



i n w hich rs(r) is given implicitly by Eqs. ( |2.2|), (^5 ), o r 
( |2.4|), w hile a(r) is given implicitly by Eqs. (2.7), ( 2~8]) , 
or ( ^.9[ ). These equations are compatible with Eq. ( 2.12 ) 
provided that the relations 



and 



M=ics 3 (xo) = ypaV(xo) 



(2.13) 



(2.14) 



hold. Equation ( p. 14 ) implies that the mass of the black 
hole is equal to the mass removed from the interior of S , 
as w as sai d previously. It can be checked that Eqs. ( 2.12 ) 
and fl2.13p ensure the continuity of the extrinsic curvature 
across E: the union of the Schwarzschild and FLRW met- 
rics forms a g lobally valid solution to the Einstein field 
equations |23|. 



B. Double- null coordinates 

The metric of the Einstein-Straus spacetime has thus 
far been expressed in terms of two coordinates patches. 
We now construct a single coordinate system (u,v,9,(j)) 
that covers both regions of the spacetime. These coor- 
dinates are null, and in the cosmological region they are 
defined by 



n-x, 



77 + x- 



(2.15) 



They are related to similar coordinates that can be in- 
troduced in the Schwarzschild region, 



t 



t + r* 



(2.16) 



where r* = J f- 1 dr = r + 2M ln(r/2M - 1). The rela- 
tions u(u), v(v) can be obtained by using these coordi- 
nates to express continuity of the induced metric across 
E. 

For example, equating the ds 2 -, of Eq. ( 2.11), expressed 
in terms of dv, to the ds 2 ^ of Eq. ( |2.10| ), expressed in 
terms of dv, yields 



(1 - 2M/r T )dv 2 - 2drxdv = a 2 dv 2 . 



(2.17) 
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Here, the scale factor must be expressed as a function of 
v — xoj an d ty, is view ed as the function of v defined by 
Eq. Q2.12D - Equation J2.17| ) is a qu adratic equation for 
dv/dv. Solving this, using Eq. (2.14), gives 



dv 
dv 



a's 



■ a a/1 — ksQ 2 



Cs 2 /a 



(2.18) 



where a prime indicates differentiation with respect to 
v ~ XOi so = s(xq ), an d C = 87rpq 3 /3. More expl icitly , 
substituting Eqs. (fD|), (|^), and ( g3| ) into Eq. ( |2.1^ ) 
yields 



7\/ 



sin xo 



[l - cos(i; - xo)] [cos xo - cos v] 
cos(2xo) - cos(w - xo) 



Q2.14D into Eq. ( |2.18[ ) yields 
dv M [cosh(u - xo 



f or k = 1. Simi larly , substituting Eqs. (|2.6| ), ( 

l] [cosh v — coshxo] 



(2.19) 
), and 



dv sinh 3 xo 



cosh(u - xo) - cosh(2xo) 



for k 



-1. Finally, using Eq. (2.9) instead gives 
M (v-xo) 3 



dv 
dv 



2xc 



3xo 



(2.20) 



(2.21) 



for k — 0. These relations for dv/dv go to zero at v — xo, 
which represents the past singularity of the cosmological 
spacetime, and they are singular at v = 3xo ; which rep- 
resents the extension of the black hole's past horizon into 
the cosmological region (see Fig. 1). 

A nearly identical set of calculations gives us the rela- 
tions u(u). Here we find 



du 
du 



M [l - cos(u + xo)] [cos xo - cos u] 



sin xo 



for k = 1, 

du 
du 



M 



cos(2xo) - cos(w + xo) 
[cosh(w + xo) — l] [coshw — coshxo 



(2.22) 



sinh xo 



for k = — 1, and 

du 
du 



cosh(u + xo) - cosh(2xo) 



M (u + xof 



2xo d u + 3xo 



(2.23) 



(2.24) 



for k = 0. 

Equations (§l|)-(f!2j) can all be integrated in closed 
form. For k = ±1, the resulting expressions for v(v) and 
u(u) are quite complicated, and we shall not display them 
here. For k = the relations are simple, and we find 



M 



41n(^-3 



H + 81n2 
3 



(2.25) 




cosmological retarded time 

FIG. 2. The relationship between the cosmological re- 
tarded-time coordinate u and the Schwarzschild retarded-time 
coordinate u. The plot is constructed for the interval 
16 < u < 30, choosing xo — 1 and setting 2M = 1. A 
graph of v(v) would be almost identical for these choices of 
parameters. 



and 

u = M 



1 / u 
6 VXo 



-41n ^ 
VXo 



(2.26) 



We have tuned the const ants of int egrat ion to produce 
agreement between Eqs. ( |2.25| ) and ( |2.26| ) and the defin- 



ing relation 



2Mln(r/2M - 1); the right- 



hand side must be evaluated on £ and expressed in terms 
of u and v by means of Eq. ( 2.12| ). The relation u(u) is 
plotted in Fig. 2; the relation between v and v is nearly 
identical. 

Both sides of the Einstein-Straus spacetime can be 
covered by the single coordinate system (u 7 v 7 9,(f>) con- 
structed above. Our coordinates, however, do no extend 
beyond the future horizon of the Schwarzschild region. 
As far as we are aware, this coordinate system was never 
presented before in the literature. 



III. SCALAR WAVE EQUATION 

We consider the propagation of a scalar field $ in the 
Einstein-Straus spacetime. This field satisfies the wave 
equation 



(□ - £R)$ = 0, 



(3.1) 



in which □ = g a/3 'V a 'Vp is the curved-spacetime 
d'Alembertian operator, R the Ricci scalar, and £ an 
arbitrary dimensionless constant; the scalar field is mini- 
mally coupled to curvature if £ = 0, and it is conformally 
coupled if £ = |. 

In th e co smological region, where the metric is given 
by Eq. (|2.5|), we separate the variables by writing 
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Mv>x)Y lm (0,4>). (3.2) 



ini 



Here, Yi m (9,<p) are the usual spherical harmonics, and 
TpiiVi X) is a reduced wave function that satisfies 



drf 



dx 2 



2a(r,)/C s 2 ( X ) 



$\ = 0, (3.3) 



where 2a(rj)jC and s(x) are respectively equal to 1 — cost? 
and sin x if k = — 1 , cosh 77 — 1 and sinh x if fc = 1 , and 



hr] 2 and x if /c = 0. 



In th e Sc hwarzschild region, where the metric is given 
by Eq. we let 



$(t,r,6,cf>) = - J2Mt,r)Y lr 



The reduced wave function now satisfies 



(3.4) 



dt 2 



dr* 2 



f 



1(1 + 1) 



2M 



fl>l = 0, (3.5) 



where r* = r + 2Mhx(r/2M - 1) and / = 1 - 2M/r. 

The field $ and its derivatives <I> iQ are continuous 
across the boundary separating the Schwarzschild and 
cosmological regions. Smoothness of the metric across X 
ensures that the reduced wave function ipi and its deriva- 
tives are also continuous at the boundary. 

Using the coordinates u and v introduced in Sec. II B, 
the reduced wave equation becomes 



dudv 



~V t (u,v) ipi 



(3.6) 



on both sides of the boundary. The effective potential 
Vi(u, v) takes a different form in the two portions of the 
spacetime. In the cosmological region, it is 



V, out = 



1(1+1) l-6£ 



s 2 ( X ) 2afa)/C" 



where 77 = ^(v + u) and x — \(v — u). 
Schwarzschild region, we have instead 



1 du dv 



2M 
r 



1(1 + 1) 



2M 



(3.7) 



In the 



(3.8) 



the 



in which r(u, v) is defined implicitly by r* = \(v— u : 
relations u(u) and v(v) were obtained at the end of Sec. II 
B. It should be noted that the potential is discontinuous 
at the boundary; this is explained by the fact that the 
density p is discontinuous at S. 

The reduced wave equation, in the form of Eq. (3.6), 



can be straightforwardly integrated on both sides of S, 
which is described by the relation v = u+2xs)- The initial 
data required for a unique solution consists of ipi(u,v = 
Vo) and tpi(u — Uq,v), the value of the wave function on 
the two null surfaces v — vq and u — uo', these values can 
be specified freely. 



IV. NUMERICAL RESULTS 

We now present our numerical integrations of Eq. . 
The algorithm used for this task is based on the finite- 
difference scheme suggested by Gundlach, Price, and 
Pullin pi . In this scheme, the (u, y ) continuum is dis- 
cretized in units of A, and Eq. ( |3.6| ) is replaced by the 
approximation 



ipi(u + A,v + A) 



1, 



1-^A 2 + 0(A 4 ) 

o 



ipi(u + A,v) 
+ ipi(u,v + A) -ipi(u,v), (4.1) 



in which the potential is evaluated at the off-grid point 
(u + iA, v + iA). For simplicity we restrict our atten- 
tion to spatially-flat cosmologies (k = 0), and a scalar 
field th at is minimally coupled to curvature (£ = 0). 
In Eq. (3/7) we must therefore substitute s(x) = X. = 
%(v-u) and 2a(rf)/C = \tf = |(« + u) 2 . In Eq. jfijh, 
the r elatio ns v(v) and u(u) are obtained from Eq. ( |2.25| ) 
and 
Eqs. 



( |2.26|) , resp ective ly; their derivatives are given by 
( |2.21| ) and ( 2.24 ). For the initial value of the wave 
function on u — uq we choose a gaussian wave packet of 
width a centered at v — v c : 



%pl(u = u , v) = exp 



(v ~ Vc) 2 
2a 2 



(4.2) 



vq we make the choice ipi(u,v = vq) = 
^ 0: the initial field is a purely ingoing wave 



On v = 

i>i(u Q ,v 

packet. 

The parameter spa ce associated with the numerical in- 
tegration of Eq. ( |3.6| ) is rather large, but moving in this 
space does not affect the qualitative aspects of the wave's 
evolution. For concreteness we set xo = 1, 2M = 1, 
uq = 17, vq — 16, v c = 16.9, a — 0.2. These choices 
ensure that the initial wave packet is contained entirely 
within the Schwarzschild region of the spacetime; its cen- 
ter is located in the strong-field region near r* — 0. The 
discretization unit A must be set sufficiently small to en- 
sure that the potential Vi in Eq. ( |4.lD is well sampled. 
In our simulations we use A = 0.001; choosing a smaller 
value does not appreciably change our results. In Fig. 3 
we display the numerical grid and a few relevant features 
of the spacet ime . 

Equation ( |4.l| ) is used to numerically evolve ipi from 
the initial conditions imposed on the lines u = 17 and 
v = 16. The wave function is then evaluated on three 
different curves of the (u, v) plane. The first is the 
line u = 30, which approximates the black hole's event 
horizon (see Fig. 3); there ipi is expressed as a func- 
tion of cosmological advanced time v. The second is 
the line v = u + 2, which represents the boundary E; 
there ipi is expressed as a function of conformal time 
rj = v — l = ii+l. The third is the line v = 30, which 
approximates future null infinity (see Fig. 3); there tpi is 
expressed as a function of cosmological retarded time u. 
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Scalar field at event horizon: monopole moment 
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FIG. 4. A plot of (u = 30)|, the absolute value of the 
I = wave function on the event horizon, as a function of 
cosmological advanced time v. Prior to v = 19, the wave be- 
haves as it would in Schwarzschild spacetime. This early-time 
behavior is characterized by quasi-normal oscillations (which 
on a logarithmic scale are revealed by deep troughs) and the 
onset of inverse power-law falloff. At v = 19 the field comes 
in causal contact with the cosmological region and suddenly 
alters its behavior. The wave function changes sign, and it 
subsequently falls off at a much slower rate; as we show in 
Sec. V, this late-time behavior is not well described by an 
inverse power-law. 



FIG. 3. The portion of the Einstein-Straus spacetime that 
is relevant for our numerical integrations. With \o = 1 5 this 
corresponds to the intervals 17 < u < 30 and 16 < v < 30. 
The thick, dashed line described by v — u + 2 represents 
the boundary £. The Schwarzschild region of the spacetime 
appears above the boundary, while the FLRW region appears 
below. The solid lines are curves r*(u, v) — constant, with the 
label giving the value of the constant (in units 2M — 1); these 
curves are meaningful only in the Schwarzschild part of the 
spacetime. Notice that the curve r* = is well approximated 
by u = u. The initial wave packet on it = 17 is centered 
at v = 16.9, at which r* ~ 0; it is therefore located in the 
strong-field region of the Schwarzschild spacetime. Because r* 
is negatively large everywhere on u = 30, this surface closely 
approximates the black hole's event horizon. We similarly 
take v — 30 to approximate future null infinity, although the 
approximation becomes poor as u increases. 



The results of our numerical simulations are displayed in 
Figs. 4-12, and the captions describe the graphs in de- 
tail. In the following paragraphs we describe and explain 
the salient features. 

Wave propagation in the Einstein-Straus spacetime is 
gove rned by the potential Vi displayed in Eqs. (3.7) and 
( |3.8[ ). The most important aspects of this potential are 
that it peaks near r* — (because of the strong cur- 
vature near the black hole) and that it is discontinuous 
at the boundary (because of the jump in the density p). 
For the purposes of understanding the gross properties 
of the wave's evolution, the potential can be thought of 
as consisting of a narrow barrier near r* = and a jump 
discontinuity at x = Xo- After each encounter with the 
barrier or the jump, the wave is partially transmitted and 
reflected, and these processes determine the late-time be- 
havior of the scalar field. The diagram of Fig. 3 informs 
us that the curve r* = is well approximated by the line 
v = u. The potential barrier is therefore located near 
v = u, while the jump discontinuity is at v = u + 2. 

Consider, as we do in our numerical simulations, an 
initially ingoing wave packet on u = 17, centered at v — 
16.9. In the absence of any potential, the wave would 
keep its shape as it moves toward increasing values of 
u. In the diagram of Fig. 3, the wave packet would move 
vertically upward, and it would eventually register on the 
"event horizon" at u = 30; the wave profile on u = 30 
would be identical to the initial configuration on u = 
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Scalar field at boundary: monopole moment 
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cosmological advanced time 

FIG. 5. A plot of \ipi(u = 30) |, the absolute value of the 
I = 1 wave function on the event horizon, as a function of cos- 
mological advanced time v. The early-time behavior is again 
characterized by quasi-normal oscillations (with a higher fre- 
quency and a faster decay time than the I = case) and the 
onset of inverse power- law falloff. At v = 19 a faint echo of 
the early-time signal registers on the event horizon. A second 
echo also appears near v = 21. Beyond v = 22, the wave 
function drops below 1CF 15 and our numerical results become 
unreliable. 
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conformal time 

FIG. 7. A plot of \ipo(v — u + 2)|, the absolute value of 
the I — wave function on the boundary E, as a function of 
conformal time rj. At early times the field behaves as it would 
in Schwarzschild spacetime. At later times, from r) = 19 on- 
ward, the multiple reflections from the potential barrier near 
r* = and the jump discontinuity at E alter the behavior of 
the wave function. At rj ~ 19.7, the wave function changes 
sign, and it subsequently falls off at a slow rate; as we show 
in Sec. V, this behavior is not well described by an inverse 
power-law. 



Scalar field at boundary: dipole moment 



Scalar field at event horizon: quadrupole moment 




cosmological advanced time 

FIG. 6. A plot of \^)i(u — 30)|, the absolute value of the 
I = 2 wave function on the event horizon, as a function of 
cosmological advanced time v. Again we have quasi-normal 
oscillations and inverse power-law falloff at early times, fol- 
lowed by echos of the early signal at later times. Our numer- 
ical results become noisy beyond v = 20, but a second echo 
can still be seen near v = 21. 
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22 

conformal time 



FIG. 8. A plot of \ipi(v = it + 2)|, the absolute value of the 
I = 1 wave function on the boundary E, as a function of con- 
formal time rj. Again we have Schwarzschild-type behavior at 
early times, and strong deviations from it starting at rj = 19. 
The multiple reflections of the wave packet give rise to the 
oscillatory behavior seen in the figure. The period of these 
oscillations is approximately Ar/ ~ 2, twice the light-travel 
time between the boundary and the potential barrier. (At 
first glance, the frequency appears to be twice as large be- 
cause the repeated pattern, which can be seen in the interval 
18 < rj < 19, already contains a sign change.) The end result 
is a field that decays faster than in Schwarzschild spacetime, 
contrary to what takes place for I = 0. Our numerical results 
become unreliable beyond rj = 24. 
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Scalar field at boundary: quadrupole moment 



Scalar field at null infinity: dipole moment 




20 21 22 23 24 25 

conformal time 

FIG. 9. A plot of \ip2(v — u + 2) |, the absolute value of 
the I = 2 wave function on the boundary E, as a function 
of conformal time rj. We recognize the late-time oscillatory 
behavior that was identified in Fig. 8. Here also the field 
decays faster than it would in Schwarzschild spacetime. Our 
numerical results become unreliable beyond rj = 24. 
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cosmological retarded time 

FIG. 11. A plot of \i[ii(v = 30)|, the absolute value of the 
I — 1 wave function at future null infinity, as a function of 
cosmological retarded time u. The behavior seen here is more 
or less a direct image of what was already seen at the bound- 
ary, in Fig. 8. Contrary to the I = case, the field decays 
faster in the Einstein-Straus spacetime than in Schwarzschild 
spacetime. 



Scalar field at null infinity: monopole moment 
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18 20 22 24 26 28 30 

cosmological retarded time 

FIG. 10. A plot of \ip (v = 30)|, the absolute value of the 
I — wave function at future null infinity, as a function of 
cosmological retarded time u. The behavior seen here is more 
or less a direct image of what was already seen at the bound- 
ary, in Fig. 7. Compared with the Schwarzschild behavior, 
the field decays at a slower rate; as we show in Sec. V, this 
behavior is not well described by an inverse power-law. 



Scalar field at null infinity: quadrupole moment 
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cosmological retarded time 



FIG. 12. A plot of \ip 2 (v = 30)|, the absolute value of the 
I — 2 wave function at future null infinity, as a function of 
cosmological retarded time u. The behavior seen here is more 
or less a direct image of what was already seen at the bound- 
ary, in Fig. 9. Compared with its Schwarzschild behavior, the 
field decays at a faster rate. 



8 



17. Instead, our wave packet moves in the presence of 
the potential, which produces partial transmissions and 
reflections as well as distortions in the packet's shape. 

Because the initial wave packet is centered near r* = 0, 
it almost immediately encounters the potential barrier as 
it propagates away from u = 17. The partially trans- 
mitted wave moves vertically upward in the diagram of 
Fig. 3, and a distorted wave profile registers on the event 
horizon; this profile still has a recognizable center at 
v = 16.9. On the other hand, the partially reflected wave 
moves horizontally toward the right, and it soon arrives 
at the boundary; the center of the wave packet arrives 
shortly after ?7 = uo + l = 18. 

At the boundary, the right-moving wave encounters 
the jump in the potential and undergoes partial trans- 
mission and reflection. The transmitted wave proceeds 
toward future null infinity without much further distor- 
tion; the center of the wave packet arrives shortly after 
u = uq = 17. On the other hand, the reflected wave 
moves vertically upward in the diagram, and the center 
of the packet follows the line v = r/ + 1 = 19. It encoun- 
ters once more the potential barrier near u = v. 

At the barrier, the up-moving wave is partially trans- 
mitted and reflected. The transmitted wave continues 
along v = 19 and it finally registers on the event horizon, 
where it is recognized as a faint echo of the original signal. 
The reflected wave moves horizontally along u = 19, and 
proceeds toward the jump at v — u + 2. There the entire 
cycle of partial transmissions and reflections repeats once 
more. 

In the figures, the behavior of the wave function in the 
Einstein-Straus spacetime is contrasted with the behav- 
ior it would have in Schwarzschild spacetime. Because 
the scalar field is at first causally disconnected from the 
cosmological region, the early-time behaviors coincide. 
As the field comes in causal contact with the bound- 
ary and the expanding universe beyond it, its behavior 
changes. In particular, the inverse power-law falloff wit- 
nessed at late times for Schwarzschild spacetime is not 
seen in the Einstein-Straus spacetime. Instead, the field's 
late-time behavior is quite complicated, and governed by 
many cycles of partial transmissions and reflections, as 
was described above. 



V. PARTIAL TRANSMISSIONS AND 
REFLECTIONS: A TOY MODEL 

In this final section we firm up the claim made in 
Sec. IV that late-time wave propagation in the Einstein- 
Straus spacetime is governed mostly by partial transmis- 
sions and reflections. We recall that these are produced 
by the most important features of the potential Vf. a 
fairly well localized barrier near r* = and a jump dis- 
continuity at the boundary E. To explore the physics 
of those transmissions and reflections, we consider a toy 
model involving a simpler potential that still captures the 



essential features of the original potential. The new wave 
equation is sufficiently simple that we will be able to find 
an analytical solution. 



A. Formulation of the toy model 



The toy model is based on the wave equation 



d 2 



tf- 



dr/ 2 d\ 2 
in which the potential is given by 



VU = 0, 





-2/ V 2 



0<X<Xo 
X > Xo 



(5.1) 



(5.2) 



The potential vanishes on the left-hand side of the bound- 
ary at x — Xqi an d the wave propagates freely in the 
interval < \ < xo- This interval is our model for the 
region between the potential barrier near r* = and the 
boundary S: we effectively "switch off" the potential in 
the Schwarzschild region of the spacetime. The partial 
reflections off the potential barrier are replaced by to- 
tal reflections at % = 0, where we impose the boundary 
condition 



^(r?,x = 0) = 0. 



(5.3) 



On the right-hand side of the boundary, for x > Xo> the 
potential is 77-dependent, a behavior that is inherited by 
taking the special case I = 0, £ = 0, and 2a( rj)/C — ^rj 2 
of Eq. (3/7). For x> Xo, therefore, Eq. ( |5.l| ) describes a 
I = 0, minimally-coupled wave propagating in a spatially- 
flat cosmology. Notice that for x < Xo, Eq. ( |5~l| ) de- 
scribes a I = wave propagating in flat spacetime. 

The toy model incorporates an infinite potential bar- 
rier at x = and a jump discontinuity at x — Xo- We 
believe that this is sufficient to capture, in the simplest 
possible way, the essential features of wave propagation 
in the Einstein-Straus spacetime. The toy model, how- 
ever, does not incorporate an event horizon, and it does 
not take into account the effect of the 1(1 + 1) part of the 
potential. 



In the rest of the section we will integrate Eq. (|5.1| ) 
analytically, subjected to the initial conditions 



ip(u, v ) = f(u), tp(u , v) = 0, 



(5.4) 



where f(u) is a specified function with support in the 
interval vq < u < uq; this corresponds to an initially 
outgoing wave packet. We have re- introduced the null 
coordinates u and i>, related to the original coordinates 
by u = i] — x, v — r\ + x- We will be interested mostly in 
the behavior of the wave function at x — Xo- We sketch 
the evolution of the wave packet in Fig. 13. 



9 




X 



FIG. 13. Wave propagation in the toy model. The diagram 
sketches the evolution of the initial wave packet, which is 
placed on the line v — vo- When it first reaches the boundary 
a t X — Xo, the wave is partially transmitted and reflected. 
The reflected wave moves back toward \ — 0, where it is 
totally reflected. The new outgoing wave encounters once 
more the boundary, and it is again partially transmitted and 
reflected. The entire cycle repeats ad infinitum. The first 
epoch of wave propagation refers to the wave's first encounter 
with the boundary; it begins at r\ — uo + xo and ends at 
r\ = U3 + 3xo- The second epoch refers to the wave's second 
encounter with the boundary; it begins at r\ = uo + 3xo and 
ends at r] = 113 + 5xo ■ 



B. Elementary solutions 

We shall impose that ip and its derivative with respect 
to x be continuous at x = Xo- This will allow us to de- 
scribe the partial transmissions and reflections that take 
place at the boundary. The t otal reflections occurring at 
X = are described by Eq. (5.3). The difficulty of the 
forthcoming calculations resides entirely with the imposi- 
tion of these matching conditions. Solving the wave equa- 
tion separately in the two domains x < Xo and x > Xo is, 
by contrast, quite trivial. For V — 0, the general solution 
to Eq. (O) is 



V>ieft = a(u) + b(v), 



(5.5) 



in which a and b are arbitrary functions of their argu- 
ments. For V = — 2/rj 2 , the general solution to Eq. (5.1) 
is 



hi. 



d 

du 



A(u) 



d 
dv 



B(v), (5.6) 



in which A and B are arbitrary functions of their argu- 
ments. 

We note in passing that Eqs. (5.5) and J5.6] ) can easily 
be generalized to potentials that include a 1(1 + l)/x 2 
term. Starting from a monopole solution ipo, a dipole (I = 
1) solution is obtained by letting ipx — (— d/d\+ l/x^o, 
and a quadrupole (/ = 2) solution is obtained by letting 



V>2 = (—d/dx + 2/x)ipi- The general rule is that a (l + l)- 
pole solution is obtained from a Z-pole solution by letting 

ip l+1 = {-d/d x + (i + i)/xWi- 



C. Matching: First epoch 



We begin by matching the solutions ( |5.5| ) and (5.6) 
in the interval u + xo < V < u o + 3xo- During this 
first epoch, the following events take place: (i) the ini- 
tially outgoing wave packet represented by f(u) arrives 
at x = Xo an d is partially transmitted and reflected; (ii) 
the reflected wave moves back toward x = and is to- 
tally reflected there; (iii) the new outgoing wave, which 
later will be represented by the function / new (u), starts to 
move toward x = Xo- The first epoch, therefore, includes 
only the first interaction between the outgoing wave and 
the boundary. 

During the first epoch, 



V'icft = f(u)+g(v), 



(5.7) 



where f(u) represents the initially outgoing wave packet, 
and g(v) its partial reflection off the boundary; we have 
g(v) = for v < uq + 2xo- On the other hand, 



bright = F(u) F(u), 

V 



(5.8) 



and this represents the wave that is partially transmitted 
across the boundary; an overdot indicates differentiation 
with respect to u. The functions g(v) and F(u) are deter- 
mined by imposing continuity of ip and tp tX across x = Xo- 
The equation that determines F(u) is 

F(u) - + -^M™ = /(„). (5.9) 

u + Xo 2(u + xor 

The equation that determines g(v) is 

g(v) = -f(v - 2 X0 ) + F(v ~ 2xo) - ^ ~ ^ - (5.10) 

v-Xo 



We must first integrate Eq. (5.S) subjected to the ini- 
tial conditions F(uq) = F(uo) = 0. Standard integration 
techniques [E4| reveal that the required solution is 



F(u) = -(u + Xo )^' 



/(«') 



+ ^(^ + Xo)^ 7 



u {U' + X0) 1+1 
/(«' 



du' 



K + xo) 1 " 7 



— du', (5.11) 



where 7 = \j Substituting this into Eq. (5.10) gives 



9(v) = ^(w-Xo) 7 / -T-, 

2 Ju +2 X0 (v '-Xo) 1+7 



7 



(v-Xo) 7 



' ^'. (5.12) 

, +2x „ K-X0) 1 - 7 
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We recall that this ingoing wave is totally reflected at x 
0, and that it becomes the new outgoing wave f n ew(u) 
-9{v). 



By substituting Eq. fl5.1l|) into Eq. (gj) or Eq. fljUg ) 
into Eq. (5.7), we finally obtain the wave function ip(u, v). 
Its expression on the boundary is 



TO, Xo) = f(v Xo) + ^ 7 f /(??/ /1 ^ X0) */ 

Z "'"o+Xo " 

" M-XO);., (513) 



7 



"o+Xo 



/1-7 



7 'l+7 



This expression is valid in the interval uq + Xo < ?7 < 
mo + 3xo only. 



Scalar field at boundary: monopole moment 
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D. Matching: Second epoch 

We continue our int egra tion o f Eq . ( |5.1| ) and carry out 
the matching of Eqs. (5J) and (plj) in the interval uq + 
3xo < r/ < uq + 5xo- During this second epoch, the 
following events take place: (i) the new outgoing wave 
fnew(u) arrives at x — Xo an d is partially transmitted 
and reflected; (ii) the reflected wave moves back toward 
X = and is totally reflected there; (iii) the new outgoing 
wave starts to move toward x — Xo- 

During the second epoch, 



^left = /ncw(u) + 9new(v), 



(5.14) 



where 



7, 



/new(u) = --(u-XoV I , , u+7 

± Ju +2 Xo («' - Xo) 1+7 



/K-2xo) 



+ |(«-Xo)" 7 



/K-2Xo).^ 



1+2xo («' - Xo) 1 " 7 



(5.15) 



represents the new outgoing wave, and g nevr (v) its partial 
reflection off the boundary. On the other hand, 



bright = F(u) - -F(u) 
V 



(5.16) 



represents the partially transmitted wave; because F(u) 
is continuous at u — uq + 2xo, there is no need to distin- 
guish the second-epoch F(u) with a label "new". 

Th e fu nctio ns F (u) and g n ew(v) are determined by 
Eqs. (5.9) and ( 5.10| ), in which we replace f(u) by / MW (u) 
and g(v) by g n cw{v)- Equation ( |5.9| ) is solved under the 
conditions that F(uq + 2x o) an d F(u + 2xo) match the 
values obtained from Eq. (5.11). This gives 



F(u) = -(u + xo) 1+ ' 1 



2 

+ 



/(«') 



K + Xo) 1+7 



du' 



/new(^ ) 
„o+2xo K + Xo) 1+7 



conformal time 

FIG. 14. A plot of the absolute value of the wave function 
on the boundary, as a function of conformal time W- The soli d 
curve is obtained from a numerical integration of Eq. 



starting with the gaussian wave packet described by Eq. (4.2). 
For this integration we set no = 17, vq — 17, v c = 17.3, 
a — 0.1, and xo = !• The dashed curve is obtained from 



the analytical results of Eqs. (5.13) and (5.15), in which we 
substitute the totally reflected gaussian wave packet for f(u). 
The two curves labelled "toy model" start to diverge at the 
end of the second epoch, when r\ = uo + 5xo = 22. The dot- 
ted curve represents a monopole wave in the Einstein-Strauss 
spacetime, which displays a similar late-time behavior. 



T^ + Xo) 1 - 7 



mo+2xo 



/(«') 



/new(^ ) 



(«' + Xo) 1 " 7 



du' 



«o+2xo 



(«' + Xo) 



1-7 



du' 



(5.17) 



and Qncw M can be obtained immediately from 
Eq. J5.1C| ); because this result is not needed for further 
calc ulatio ns, we sha ll no t display it here. Substituting 
Eq. ( 5.17 ) into Eq. ( 5.16 ) and evaluating on the bound- 
ary, we obtain 



TO, Xo) = /ncw(?7 ^ Xo) + 2 7 7 7 



"0+3X0 f(r]l _ Xo) 



"o+Xo 

V /ncw(V - Xo) 

"o+3xo ^ 



Vl+7 



/1+7 



drj' 



7 -7 
2 ' 



"0+3X0 _ Xo) ^ / 



"o+Xo *7 

" /newtV - Xo) 

"o+3xo ^ 



/l— 7 



(5.18) 



This expression is valid in the interval uq + 3xo < V < 
mo + 5xo only. 
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E. Wave propagation in the toy model 



Equation fl5.13| ) gives V'lZLXo) i n the interval uq + xo < 
rj < uo + 3xo, and Eq. ( 5.18| ) does the same in the in- 
terval uq + 3%o < rj < uq + 5xo- It can be checked that 
these values are cont inuous at r\ = uq + 3xo- While in- 
tegration of Eq. ( 5T ) could easily be carried out beyond 
the second epoch by repeated application of the general 
method, we shall stop here: to continue beyond this point 
would give rise to messy expressions that wou ld not teach 
us anything new. The message of Eqs. ( 5.1 3| ) and ( 5.1§| ) 
is already clear: The multiple transmissions and reflec- 
tions that take place at x = Xo produce a wave function 
with a fairly complicated mathematical description, and 
its late-time behavior is not well described by a simple 
inverse power-law. 

In Fig. 14 we compare the wave function obtained by 
integrating Eq. ( |5.lD to a I = wave in the Einstein- 
Straus spacetime. The figure reveals that both waves 
have a similar late-time behavior, providing support for 
our claim that wave propagation in the Einstein-Straus 
spacetime is governed mostly by partial transmissions 
and reflections. 
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